Cost of material or information flow in complex transportation networks 
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To analyze the transport of information or material from a source to every node of a network we 
use two quantities introduced in the study of river networks: the cost and the flow. For a network 
with K nodes and M levels, we show that an upper bound to the global cost is Co,,„aa; oc KM . 
From numerical simulations for spanning tree networks with scale-free topology and with 10^ up 
to 10^ nodes, it is found, for large K, that the average number of levels and the global cost are 
given by M oc In(isr) and Co oc K\n{K), respectively. These results agree very well with the ones 
obtained from a mean-field approach. If the network is characterized by a degree distribution of 
connectivity P{k) oc we also find that the transport efficiency increases as long as 7 decreases 
and that spanning tree networks with scale-free topology are more optimized to transfer information 
or material than random networks. 

PACS numbers: 05.10.-a, 05-45.-a, 87.18. Sn 



There is a great interest in the study of networks and 
its applications pj, i2|, i3|| because complex network struc- 
ture underlies many biological, social and technological 
systems. An interesting development was the discov- 
ery that for many real networks, such as Internet 
metabolic networks [H] and citation networks Q, the de- 
gree distribution follows a long-tailed power law relation- 
ship P{k) oc fc~^, where k is the number of links of a node 
[tJ. On the other hand, the degree distribution of a ran- 
dom graph is given by a Poisson distribution with a peak 
at P{{k)). Note that in a random graph the edges are 
placed randomly, implying that the greater part of the 
nodes have approximately the same degree (fc) as in the 
classical Erdos-Renyi model. 

The majority of papers about complex networks are 
focused on static topological properties or on models for 
their growth. In this paper we are interested in the infor- 
mation, material or energy flow in complex transporta- 
tion networks. Using tools developed to study trans- 
portation networks @, H, , we study the relation be- 
tween network topology and efficiency to transfer infor- 
mation, material or energy from a source to every node 
of the network. Us ing a simple argument developed for 
us and a colleague [Tl| , we first determine that an upper 
bound for the global cost for a network with K nodes and 
M levels is C()_max oa MK. Then we show for a spanning 
tree network with scale-free topology that the number of 
levels grows with the number of nodes as M oc h\{K) 
and that the cost is proportional to K\x\{K). Both re- 
sults, which are valid in the limit of large if, are obtained 
from numerical simulations for large networks {K vary- 
ing from 10^ up to 10^ nodes) and from a mean-field 
approach. This implies that Co oc Co,max MK for a 
spanning tree network with scale-free topology. 

In general, the analysis of network topology is done by 



using concepts such as small world character fl2l] , cluster- 
ing coefficient [l^l and degree distribution |7|]. However, 
transport properties have not been extensively studied. 
Lopez et al. [l3| have investigated the transport prop- 
erties of scale-free and Erdos-Renyi networks by analyz- 
ing the conductance G between two arbitrarily chosen 
nodes. They showed that for scale-free networks with 
7 > 2 the conductance display a power law tail distribu- 
tion (j>sF{G) oc G~^°, that is related to the degree distri- 
bution P{k) by go = 27 — 1. On the other hand, the con- 
ductivity distribution of Erdos-Renyi networks 4>er{G) 
decays exponentially. The authors concluded that trans- 
port in scale-free networks is better than in Erdos-Renyi 
ones. Castro e Silva et al. HI] have studied a deter- 
ministic Boolean dynamics on scale-free networks using 
a damage spreading technique. They showed that the 
Hamming distance and the number of I's exhibit power 
law behavior, with the exponents depending on the value 
of 7. 

Since the optimization of information (material) flow in 
networks is an issue of great importance in diverse disci- 
plines [15], it is natural to ask which network topology is 
more efficient to transfer information (or material) from 
a source to every node of the network. To answer this 
question, we use tools developed to study transportation 
networks [1, [^, ■ A transportation system is composed 
by a source (a central node) and a set of K nodes to be 
reached. Each of the nodes is connected to one or more of 
its neighbors nodes in such way that there is a route from 
the source to every node of the network. This character- 
izes a spanning network that can present loops. Since 
any transportation network has, in general, a main route 
from the source to every node, we can define a spanning 
tree of the network as a loop less subset of the network in 
which each node can be reached from the source by the 
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flow H, [1, [T3|- Here, let us consider only spanning tree 
networks, in which the flow rate does not change with 
time and that the source has an outward flow whose ex- 
actly equals the sum of the flow rate into all the links. 
Examples of such networks include the vascular system of 
mammals [l6l |. electric circuits and river drainage basins 
i- 

To generate the spanning tree networks with scale-free 
topology, we use a growing method with re-direction of 
links [131 • We start with two nodes, the source X and a 
node Y . The node X or the source is its own ancestor 
and the ancestor of node Y . To insert a new node Z , we 
first select a node of the network with equal probability. 
With probability 1 — r the link is established between the 
node Z and the selected node, and with probability r, the 
link is established with the ancestor of the selected node. 
When we have r = (no re-direction) , a new node Z can 
be directly connected to each node of the network with 
equal probability. This constitutes a random network 
with Poisson degree distribution of connectivity. For val- 
ues of r in the interval [0, 1] , the network has a long-tailed 
scale-free distribution of connectivity with different val- 
ues of 7. Since the relation between r and 7 is given by 
7 = 1 + the case r = 0.5 corresponds to the linear 
preferential attachment (7 = 3) 0]. If r = 1, the link is 
always established with the ancestor of the selected node 
and the network has a "star-like" configuration, in which 
all the nodes are connected directly to the source. In 
this case, the source has K links and the others K nodes 
have only one link. The number of levels M of a network 
can be defined as the number of links among the source 
and the most distant node. Then, the "star-like" net- 
work has only one level, M — 1. Note that the networks 
built with this method have a tree-like structure, without 
links between nodes of the same level. Thus, each node 
is connected to the source by only one path. 

Let us measure the cost of transport in these networks 
using two quantities of the river network theory ^ and 
food- web studies [10]. The first one is the quantity Ai of 
nodes constituting the sub-tree that begins in the node 
i, plus the node i itself. The other quantity is the trans- 
portation cost, defined by the sum Ci = A^, where k 
runs over the sub tree that begins in the node i. In anal- 
ogy with river networks, Ai can be related to the water 
flux arriving in site i and Ci is identified as the volume 
of water contained in the sub tree that begins in i. La- 
behng the source by i = 0, we have that ^0 = K + 1, 
K nodes plus the source, is the total size of the network 
and Co is the global cost of the transportation system. 
For instance, let us consider two networks arrangement, 
the "star-like" and the "chain-like" configurations. In 
the first case the source is at the center and all the nodes 
are directly connected to it, while in the second case, the 
"chain-like" networks, all nodes have only one incoming 
link and one outgoing link, except the source and the 
most distant node. It is knowledge that these configura- 
tions are the extreme cases in sense of cost of transport 
from the source to every node in the network. If there is 



no constraint on the topology of the network, the span- 
ning tree can be "chain-like" , "star-like" or something in 
between. Note that the "star-like" configuration is the 
most efficient case and the "chain-like" is the less effi- 
cient one. Any other tree like structure is between these 
limits. Motivated by the works of Banavar et al. [9| and 
West et al. [l^, Garlaschelli 10] supposed that in gen- 
eral Co and K are related by a power law relationship 
Co oc if, where rj quantifies the degree of optimization 
of the transportation system. It is easy to show that for 
the "star-like" and 'chain-like" configurations we have 
that 77 = 1 and 77 = 2, respectively. 

In a previous paper, we and a colleague [ll| derived 
general relations for the upper and the lower bounds for 
the global cost Co in a general spanning tree network. 
Our argument is based in the special characteristic of 
the cost function: a node put as near as possible of the 
source has a minimum contribution to the global cost; 
and a node put as far as possible of the source has a 
maximum contribution to the global cost. So, we can un- 
derstand intuitively why the "star-like" and "chain-like" 
configurations have the minimum and maximum value of 
the cost. Let us consider a spanning tree network with M 
levels and K nodes. To obey the constraint of M levels, 
we put one node in each level. Now we must distributed 
the reminder K — M nodes. If we put all the reminder 
nodes directly linked in the source we have the network 
with the minimum value of the cost of transport. In other 
hand, if wc put the reminder nodes in the last level we 
have the network with the maximum transportation cost. 
Therefore the minimum and the maximum cost are given 
by 

M 

Co_ = 1 + 2X+— (M-1) , 

M 

Co,„.. - l+X(Af + 1) + — (1-M) . 

In the limit that K >> M and M >> 1, the above 
equations reduces to 

Co„.„ ocK , and Co„„,. oc KM . (1) 

Let us look at the two extreme configurations again. 
The mean number of levels in the "star-like" and "chain- 
like" networks scales as Mgi oc K° and Md oc K^, respec- 
tively. In the first case, using the above equations, the 
upper and lower bounds of the exponent 77 have the same 
value rjmax — Vmin = 1. lu the "chain- like" case, rj seems 
to be bounded between 1 e 2. Note that i]min = 1 is 
valid for all networks because the minimum cost is fixed. 
On the other hand, the maximum cost can vary if M 
depends on K. For a network with a topology between 
the two extreme cases, we expect that M depends on K 
and that the cost be proportional to the maximum cost 

(Co oc C'o.jnax OC KM). 

Now, let us return to the case of scale-free networks 
grown by the re-direction of links [13] . In figure ([T]) we 
have plotted the mean number of levels M against the 
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network size K for different values of r. Observe that in 
all cases we have that M oc In K. This simulation results 
suggest that in these networks the cost function is given 
by 
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FIG. 1; " (Color online)" Linear-Log plots of the mean number 
of levels versus the network size K for spanning tree networks 
with scale-free topology. The networks are built with different 
re-direction probability r. 



Co cx K\nK . (2) 

This result is in contrast with the premise that the 
shape of Cq as a function of K always follows a power-law 
relation Co oc if, where the scaling exponent 77 quanti- 
fies the degree of optimization of the transportation sys- 
tem used in the study of food webs [13, El, E, [13] ■ It 
is worth mentioning that in this problem, the number of 
levels M does not grow with K and, usually is assumed 
constant [M < 6). This implies that Co (x K {rj = 1) ioi 
food webs because the upper and lower bounds have the 
same scaling behavior. 

In order to verify if Eq. ([2]) holds for the scale-free 
topology, we evaluate numerically the global cost Co for 
networks with the size K varying from 10^ up to 10^. 
For each network size with a fixed value of the probabil- 
ity of re-direction r, we compute the mean value of the 
global cost Co{K) averaged in 10'' samples. In figure [2] 
it is shown the Co versus K plot for different values of r. 
Note that cost of transport decreases as long as r grows. 
This suggests that among the scale-free transportation 
networks the one that has the smallest value of 7 is the 
more efficient. Note that the network constructed with 
the re-direction probability r = (random network) has 
the biggest cost. Observe also that the network with 
the smallest cost of transport has the smallest number of 
levels M. Of course, for networks with the same size K, 



the ones with small number of levels M are closer to the 
"star-like" configuration than to the "chain-like" one. In 
fact, the results shown in Fig. [5] suggest that in the limit 
r 1 the slope of the line approaches 0, implying that 
Co oc K. In this case the network has a "star-like" shape 
because we have only a hub with all nodes connected to 
it. 




FIG. 2: "(Color online)" Linear-log Plots of the Co/K versus 
In K for spanning tree networks with scale- free topology. The 
networks are built with different re-direction probability r. 
The cost of transport decreases as long as r grows. 

The numerical results can be explained by a mean-field 
like approximation. Let us evaluate the average number 
of levels, A/, and the cost, Co, for a scaling-free tree 
(without loops) with K nodes. Since (k) is the average 
connectivity, then the source (node at level 0), is con- 
nected, in average, to (fc) nodes. Each one of the nodes 
directly connected to the source (nodes at level 1), are, 
in average connected to (fc) other nodes (nodes at level 
2). This argument is repeated until the level N, the last 
level, in such way that the average number of levels M is 
given by 

(fc)l + (fc)^2 + ... + (fc)^iV 
l + (fc) + (fc)2... + (fc)^ ■ 

If we define Sn = '^n=oi^)^ ' ^^"^ above equation can be 
written as M = Since the geometric series can 

be easily evaluated, we obtain that 

^ (jV+l)(fc)^+i _ (fc) 

- 1 (fc) - 1 ' 

implying that for large N and (fc) > 1 that M oc N. On 
the other hand, Sn must be equal to the number of nodes 
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K , namely 

(fc)^+i - 1 
Sn = =K . 

For N large, we have that (/c)^ oc K, implying that 

In if 



N cx 



In(fc) 



(3) 



Therefore the average number of levels behaves as M oc 
lnif/ln(fc), and agrees very well with the numerical re- 
sults (see Fig. [1]). We know that Ai is equal to the num- 
ber of nodes connected to node i plus itself. Then, the re- 
currence equation is given by Apj^n — 1 + {k)Aj^_(^n-i) , 
with the initial condition An = I. This equation can be 
solved and we obtain the solution, namely 



An-u = YM - 



i=Q 



(fc)"+l - 1 



(4) 



The transportation cost is defined by the sum Ci — 
'^f, Ak, where k runs over the sub tree that begins in the 
node i. It follows that 



C 



A 





E 

i—n — l 



Using Eq. Q, the above expression can be easily evalu- 
ated. We find that 



C 



N-n — 



(fc)"+i - 1 n(fc)"+i (1 - 



(fc)-i (fc)-i m-ir ■ 

Putting n ^ N we obtain the expression for Co, namely 



Co 



(fc)-l 



m-iy 



When N is large and (k) > 1, the above equation fur- 
nishes that Co oc {N{k}^)/{(k) - 1) . Since (fc)^ oc K 
and N oc M, it follows that Co oc MK, meaning that Cq 
is proportional to its maximum value, as supposed before 
(see Eq. ([T])). Using again Eq. ([3]), Co can be written as 



Co a 



KlnK 



((fc)-l)ln(fc) 



Moreover, this equation can be expressed as a function 
of the parameter r by using that (fc) = (7 — l)/(7 — 2) = 
1/(1 -r). The result 



Co oc -[(1 - r)is:in A']/[rln(l - r)] 



agrees very well with the numerical one (see Fig. [2]). 

In summary, we studied the transportation properties 
of random spanning tree networks and ones with scale- 
free topology. We verified that among the scale-free net- 
works, the efficiency of the transportation properties is 
related with the value of the exponent 7. Compact net- 
works, with small value of the 7, are more efficient that 
non-compact networks, with large value of 7. More- 
over, we verified also that any spanning tree network 
with scale-free topology is more efficient than random 
ones, although the global cost of both kinds of networks 
have a similar dependency (Cq oc Khi{K)) in the net- 
work number of nodes K. LAB thanks FAFESF and 
CAFES, Brazilian agencies, for partial financial support. 
JKL thanks CNFq and FAFEMIG for financial support. 
* Electronic address: lauro@ufrb.edu.br; t Electronic ad- 
dress: jaff@fisica.ufmg.br. 
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